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Abstract 

We use the fractional integrals to describe fractal media. We consider the fractal 
media as special ("fractional") continuous media. We discuss the possible experimental 
testing of the continuous medium model for fractal media that is suggested in Phys. 
Lett. A. 336 (2005) 167-174. Lhis test is connected with measure of period of the 
Maxwell pendulum with fractal medium cylinder. 
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1 Introduction 

The real structure of the media are characterized by an extremely complex and irregular ge- 
ometry [T] . Because the methods of Euclidean geometry, which ordinarily deals with regular 
sets, are purely suited for describing objects such as in nature, stochastic models are taken 
into account [2|. Another possible way of describing a complex structure of the media is to 
use fractal theory of sets of fractional dimensionality PQ. Although the fractal dimensionality 
does not reflect completely the geometric properties of the fractal, it nevertheless permits a 
number of important conclusions about the behavior of fractal structures [3 II]- For example, 
if it is assumed that matter with a constant density is distributed over the fractal, then the 
mass of the fractal enclosed in a volume of characteristic dimension R satisfies the scaling 
law M(R) ~ R D , whereas for a regular n-dimensional Euclidean object M(R) ~ R n . Let us 
assume that a medium can be treated on a scale R as a stochastic fractal of dimensionality 
D < 3 embedded in a Euclidean space of dimensionality n = 3. Naturally, in real objects 
the fractal structure cannot be observed on all scales. For example, Katz and Thompson 
jH] presented experimental evidence indicating that the pore spaces of a set of sandstone 
samples are fractals in length extending from 10 angstrom to 100 [im. 

In the general case, the fractal media cannot be considered as continuous media. There 
are points and domains that are not filled of the medium particles. We suggest to consider the 
fractal media as special (fractional) continuous media. We use the procedure of replacement 
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of the fractal medium with fractal mass dimension by some continuous medium that is 
described by fractional integrals. This procedure is a fractional generalization of Christensen 
approach 6j. Suggested procedure leads to the fractional integration and differentiation to 
describe fractal media. The fractional integrals allow us to take into account the fractality of 
the media. In order to describe the fractal media by continuous medium model, we must use 
the fractional integrals. The order of fractional integral is equal to the fractal mass dimension 
of the media. More consistent approach to describe the fractal media is connected with the 
mathematical definition of the integrals on fractals. In Ref. jEj was proved that integrals 
on net of fractals can be approximated by fractional integrals. Therefore, we can consider 
the fractional integrals as approximation of the integrals on fractals [7j. In Ref. [H], we 
proved that fractional integrals can be considered as integrals over the space with fractional 
dimension up to numerical factor. To prove this statement, we use the well-known formulas 
of dimensional regularizations [TU] . 

We use the fractional integrals to describe fractal media. We consider the simple example 
of application of fractal integrals to describe the fractal medium. We discuss the possible 
experimental testing of the continuous medium model for fractal media that is suggested in 
[T2] . This test is connected with measure of period of the Maxwell pendulum with fractal 
medium cylinder. 

It is known that the homogeneous medium cylinder with mass M and the radius R such 
that the medium of cylinder has the integer mass dimension, has the moments of inertia 

lf ] = (1/2)MR 2 , 

where z is a cylinder axis. If the fractal medium cylinder has the mass M, radius R, and 
the fractal mass dimension, then the moment of inertia has the form 

= ^—MR\ (1) 
a + 2 

where a is a fractal mass dimension of cross-section of the cylinder (1 < a < 2). The 
parameter a can be calculated by box counting method for the cross-section of the cylinder. 

As the result, we get that two homogeneous cylinders with the same masses and ra- 
diuses have the unequal moments of inertia. These cylinders have the different properties as 
Maxwell pendulums. The periods of oscillation for these pendulums are connected by the 
equation 

( Tt\2 4(q + l) 
V Tq ( 2 W 3(a + 2)" 

where is a period of oscillation of fractal medium Maxwell pendulum. For example, the 
experiment can be realized by using the sandstone cylinder. Note that Katz and Thompson 
[5] presented experimental evidence indicating that the pore spaces of a set of sandstone 
samples are fractals and are self-similar over three to four orders of magnitude in length 
extending from 10 angstrom to 100 jt/m. 
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2 Fractal Media and Fractional Integrals 



The cornerstone of fractals is the meaning of dimension, specifically the fractal dimension. 
Fractal dimension can be best calculated by box counting method which means drawing a 
box of size R and counting the mass inside. The mass fractal dimension can be easy 

measured for fractal media. The properties of the fractal media like mass obeys a power law 
relation 

M ~ R D , (D < 3), (2) 

where M is the mass of fractal medium, R is a box size (or a sphere radius), and D is a 
mass fractal dimension. Amount of mass of a medium inside a box of size R has a power 
law relation (J2J). 

The power law relation (J2J) can be naturally derived by using the fractional integral. In 
Ref [12] . we prove that the mass fractal dimension is connected with the order of fractional 
integrals. Therefore the fractional integrals can be used to describe fractal media with non- 
integer mass dimensions. 

Let us consider the region Wa in 3-dimensional Euclidean space E 3 , where A is the 
midpoint of this region. The volume of the region Wa is denoted by V(Wjl). If the region 
Wa is a ball with the radius Ra, then the midpoint A is a center of the ball, and the volume 
y(^) = (4/3)tcR\ . The mass of the region Wa in the fractal medium is denoted by 
M(Wa)- The fractality of medium means that the mass of this medium in any region Wa of 
Euclidean space E 3 increase more slowly that the volume of this region. For the ball region 
of the fractal medium, this property can be described by the power law (J2J), where R is the 
radius of the ball Wa- 

Fractal media are called homogeneous fractal media if the power law (J2J) does not depend 
on the translation of the region. The homogeneity property of the medium can be formulated 
in the form: For all regions Wa and Wb of the homogeneous fractal medium such that the 
volumes are equal V(Wa) = V(Wb), we have that the masses of these regions are equal 
M(Wa) = M(Wb)- Note that the wide class of the fractal media satisfies the homogeneous 
property. In many cases, we can consider the porous media OUS], polymers [THI, colloid 
aggregates [T7j, and aerogels [TH] as homogeneous fractal media. 

In Ref. [12j . the continuous medium model for the fractal media was suggested. Note 
that the fractality and homogeneity properties can be realized in the following forms: 

(1) Fractality: The mass of the ball region W of fractal medium obeys a power law relation 

M D (W) = M A D , (3) 

where D < 3, R is the radius of the ball, and R c is the characteristic value. In the general 
case, we have the scaling law relation M D (XW) = \ D M D (W), where XW = {Ax, x G W}. 

(2) Homogeneity: The local density of homogeneous fractal media is translation invariant 
value that have the form p(x) = po = const. 
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These requirements can be realized by the fractional generalization of the equation 



M 3 (W) = / p(x)dV 3 . (4) 



w 



Let us define the fractional integral in Euclidean space E 3 in the Riesz form The 
fractional generalization of Eq. (J3J) can be realized in the following form 



M D (W) = / P (x)dV D , (5) 



w 



where dVo = C3(-D, x)<iV3, and 



c 3 ^,x) r(B/2) |x| . 

Here we use the initial points in the fractional integrals are set to zero. The numerical factor 
in Eq. (|SJ) has this form in order to derive usual integral in the limit D — > (3 — 0). Note that 
the usual numerical factor ^(D) = T(l/2)/2 D ir 3/2 T(D/2), which is used in Ref. 02] leads 
to 7 3 ~ 1 (3 - 0) = r(l/2)/2 3 vr 3 / 2 r(3/2) in the limit D -> (3 - 0). 

In order to have the usual dimensions of the physical values, we can use vector x, and 
coordinates dimensionless values. For example, we can define x = Xi/x , y = x 2 /x , 

z = xz/xo, where xi, x 2 , and X3 are the coordinates with the usual physical dimension. 

We can rewrite Eq. (JHJ) in the form 

M D {W) = 23 "^ (3 ( 2) / p(x)|x| D - 3 rf 3 x. (6) 



T(D/2) Jw 

If we consider the homogeneous fractal medium (p(x) = po = const) and the ball region 
W = {x : |x| < R}, then we have 

Using the spherical coordinates, we get 



T(D/2) ^Jw 1 1 11 DT(D/2) 

As the result, we have M(W) ~ R D , i.e., we derive Eq. (JSJ) up to the numerical factor. 
Therefore the fractal medium with non-integer mass dimension D can be described by frac- 
tional integral of order D. 

Note that the interpretation of the fractional integration is connected with fractional 
dimension This interpretation follows from the well-known formulas for dimensional 
regularizations [TO] : The fractional integral can be considered as a integral in the fractional 
dimension space up to the numerical factor T(D/2)/ \2tt d / 2 T(D)). 
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3 Moment of Inertia of Fractal Medium Cylinder 

In this section, we consider the fractional generalization of the equation for moment of inertia. 

The equation for the moment of inertia of homogeneous cylinder with integer mass di- 
mension has the well-known form 

4 2) = Pof (x 2 + y 2 )dS 2 [ dz. (7) 

JS JL 

Here z is the cylinder axis, and dS 2 = dxdy. The fractional generalization of Eq. (J7J) for 
moment of inertia can be defined by the equation 

1^ = po Jy + y 2 )dS a J l dl p , (8) 

where we use the following notations 

dS a = c(a) \Jx 2 + y 2 ) dS 2 , dS 2 = dxdy, c(a) = , dip = dz. (9) 

V v / U«/ 2 J HP) 

The numerical factor in Eq. (jHJ) has this form in order to derive usual integral in the limit 
a — » (2 — 0) and /3 — > (1 — 0). The parameters a and (3 are 

1 <a<2, 0</3<l. 

If a = 2 and (3 = 1, then Eq. (jHJ) has form ((7j). The parameter a is a fractal mass dimension 
of the cross-section of cylinder. This parameter can be easy calculated from the experimental 
data. It can be calculated by box counting method for the cross-section of the cylinder. 
Substituting Eq. © in Eq. (jHJ, we get 



r(/3) 

Here we consider the cylindrical region W that is defined by the relations 

L = {z : < z < H}, S = {(x,y) : < x 2 + y 2 < R 2 }. (10) 
Using the cylindrical coordinates (<f),r,z), we have 

dS 2 = dxdy = rdrd<f), {x 2 + y 2 ) a/2 = r a . 
Therefore the moment of inertia is defined by 

As the result, we have 

j{a) = 27rp c(q) R g+2 H 



If a = 2 and = 1, we get I® = (l/2)irp R i H . 

The mass of the usual medium cylinder (fTUJ) is defined by the equation 

M = po [ dS 2 I dz = 2np [* rdr [ dz = np R 2 H. (12) 
Js Jl Jo Jo 

We can consider the fractional generalization of this equation. The mass of fractal medium 
cylinder (fTUj) can be defined by the equation 

M a = p J s dS Q J L dlp. (13) 

Using the cylinder coordinates, we get 

2 W (q) /V-Mr l H z^dz = 27rPoC(a) ^^ 
T(P) Jo Jo aPT(P) U " 

As the result, we have 

Substituting the mass (|14j) in the moment of inertia (jllj) . we get the relation 

JW = < -——M a R 2 . (15) 
a + 2 

Note that Eq. (fTo^) does not have the parameter f3. If a = 2, we have the well-known relation 
1^' = {1/2)MR 2 for the homogeneous cylinder that has the integer mass dimension D = 3 
and a = 2. 

If we consider the fractal medium cylinder with the mass and radius that are equal to 
mass and radius of the homogeneous medium cylinder with integer mass dimension, then 
the moments of inertia of these cylinders are connected by the equation 

jW = J^L/fO. (16) 

a+2 z v ; 

Here 1^ is the moment of inertia for the cylinder with integer mass dimension D = 3 and 
a = 2. For example, the parameter a = 1.5 leads us to the relation f( 3 / 2 ) = (6/7)1^ . Using 
1 < a < 2, we have the general relation 

(2/3) < /(«74 2) < I- 



4 Equation of Motion for Maxwell Pendulum 

Maxwell pendulum is used to demonstrate transformations between gravitational potential 
energy and rotational kinetic energy. Wind the string up on the small axis, giving the device 
some initial gravitational potential energy. When released, this gravitational potential energy 
is converted into rotational kinetic energy, with a lesser amount of translational kinetic 
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energy. We consider the Maxwell pendulum as a cylinder that is suspended by string. The 
string is wound on the cylinder. 

The equations of motion for Maxwell pendulum have the form 

M a ^ = M a g-T, lW*£ = KT, 
dt y 2 dt 

where g is the acceleration such that g ~ 9.81(m/s 2 ); the axis z is a cylinder axis, T is a 
string tension. Using v y = u z R, we have 

„ r dv v , , li a > dv v 
M—^ = Mg- ' 7 y . 
dt y R 2 dt 

As the result, we get the acceleration of the cylinder 

= (k)y_ = Mg = g 
ay dt M + li a) /R 2 1 + Iz^ /(MR 2 ) ' 

Substituting Eq. (|T5j) in Eq. (|17|). we get 

« + 2 

a„ = o. ( 18) 

For the fractal mass dimension of cross-section of cylinder a = 1.5, we get a y = (3/5)g ~ 
6.87 (m/s 2 ). For the cylinder with integer mass dimension of the cross-section (a = 2), we 
have a y = (2/3)g ~ 6.54 (m/s 2 ). The period T of oscillation for this Maxwell pendulum is 
defined by the equation 



T = 4t = 4 A /2L/a y , 

where L is a string length, and the time to satisfies the equation a y t\/2 = L. Therefore, we 
get the relation for the periods 



rp(a) . _ 

io) = V a V I a V 



4(a + l) 



lf> ~ v " ' " _ ^3(a + 2)- 
Using 1 < a < 2, we can see that 

(8/9) < (T W /T (2) ) 2 < 1. 

Note the parameter a can be calculated by box counting method for the cross-section of the 
cylinder. For a = 1.5, we have (Tq^/Tq^ = 0.952. 
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5 Conclusion 



In the general case, the fractal media cannot be considered as continuous media. There are 
points and domains that are not filled of the medium particles. In Ref. [O] . we suggest to 
consider the fractal media as special ("fractional") continuous media. We use the procedure 
of replacement of the medium with fractal mass dimension by some continuous medium 
that is described by fractional integrals. This procedure is a fractional generalization of 
Christensen approach Suggested procedure leads to the fractional integration and differ- 
entiation to describe fractal media. The fractional integrals are considered as approximation 
of the integrals on fractals [S|. Note that fractional integrals can be considered as integrals 
over the space with fractional dimension up to numerical factor 12II1 122] • To prove we use 
the well-known formulas of dimensional regularizations [TO]- The fractional integrals allow 
us to take into account the fractality of the media. 

The fractional continuous models of fractal media can have a wide application. This 
is due in part to the relatively small numbers of parameters that define a random fractal 
medium of great complexity and rich structure. In many cases, the real fractal structure of 
matter can be disregarded and the medium can be replaced by some "fractional" continu- 
ous mathematical model. In order to describe the media with non-integer mass dimension, 
we must use the fractional calculus. Smoothing of the microscopic characteristics over the 
physically infinitesimal volume transform the initial fractal medium into "fractional" con- 
tinuous model that uses the fractional integrals. The order of fractional integral is equal to 
the fractal mass dimension of the medium. The fractional continuous model allows us to 
describe dynamics for wide class fractal media [T3J I2D I2HI • 

In this Letter we consider the simple experiment that allows us to prove the fractional 
continuous media model |T2| for fractal media. This simple experiment prove that the 
fractional integrals can be used to describe fractal media. 

The fractal medium cylinder with the mass M, radius R, and the fractal mass dimension, 
has the moment of inertia in the form 



where a is a fractal mass dimension of cross-section of the cylinder (1 < a < 2). The 
parameter a can be calculated by box counting method for the cross-section of the cylinder. 
As the result, we get that two homogeneous cylinders with equal masses and radiuses can have 
the unequal moments of inertia. These cylinders have the different properties as Maxwell 
pendulums. The periods of oscillation for these pendulums are connected by the equation 



where T is a period of oscillation of fractal medium Maxwell pendulum. For example, 
the experiment can be realized by using the sandstone. Note that Katz and Thompson [3] 
presented experimental evidence indicating that the pore spaces of a set of sandstone samples 
are fractals and are self-similar over three to four orders of magnitude in length extending 



/■(«) 



a 



MR 2 , 



a + 2 
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from 10 angstrom to 100 \im. The deviation T from T is no more that 6 per cent. 
Therefore the precision of the experiments must be high. 
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